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ABSTRACT 

We present a study of the morphology of the ejecta in Supernova 1987A based on images 
and spectra from the HST as well as integral field spectroscopy from VLT/SINFONI. The HST 
observations were obtained between 1994 - 2011 and primarily probe the outer hydrogen-rich 
zones of the ejecta. The SINFONI observations were obtained in 2005 and 2011 and instead 
probe the [Si I]/[Fe II] emission from the inner regions. We find a strong temporal evolution 
of the morphology in the HST images, from a roughly elliptical shape before ~ 5, 000 days, 
to a more irregular, edge-brightened morphology with a 'hole' in the middle thereafter. We 
demonstrate that this transition is a natural consequence of the change in the dominant energy 
source powering the ejecta, from radioactive decay before ~ 5, 000 days to X-ray input from the 
circumstellar interaction thereafter. The [Si I]/[Fe II] images display a more uniform morphology, 
which may be due to a remaining significant contribution from radioactivity in the inner ejecta 
and the higher abundance of these elements in the core. Both the Ha and the [Si I]/[Fe II] line 
profiles show that the ejecta are distributed fairly close to the plane of the inner circumstellar 
ring, which is assumed to define the rotational axis of the progenitor star. The Ha emission 
extends to higher velocities than [Si I]/[Fe II], as expected from theoretical models. There is no 
clear symmetry axis for all the emission and we are unable to model the ejecta distribution with 
a simple ellipsoid model with a uniform distribution of dust. Instead, we find that the emission 
is concentrated to clumps and that the emission is distributed somewhat closer to the ring in 
the north than in the south. This north-south asymmetry may be partially explained by dust 
absorption. We compare our results with explosion models and find some qualitative agreement, 
but note that the observations show a higher degree of large-scale asymmetry. 
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1. Introduction 

Both observations and numerical simulations 
show that supernova (SN) explosions are asym- 
metric. In particular, simulations show that large- 
scale instabilities are likely to p lay a key role 
in tr i ggering the exp l osions (e.g. Kifonidis et al 

20121) 



2006t iHammer et all I2OIOI: iMiiller et a" 



The instabilities are expected to be reflected in 
the morphology of the ejecta in the homologous 
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phase, which thus provides one of the few direct 
diagnostics of the explosion mechanism. 

Information about the morphology of core- 
collapse SNe can be obtained in several different 
ways. Direct imaging of Galactic SN remnants can 
give detailed information about both the abun- 
dance and the density distributions. However, 
there are very few cases unaffected by the in- 
teraction with the circumstellar and interstellar 
media, or a central pulsar wind nebulae. By far, 
the most important case is Gas A, where both 



torial ring and the interact ion caused the first so- 



X-ra y fe.g. .Hwang fc Lamind 2012 ) and optical 
(e.g. Fesen et al.l 20061 ) observations have pro- 



vided important information about instabilities 
and large-scale mixing in the explosion. In addi- 
tion, light echoes of Gas A have revealed evidence 
for asymmetries and shown that the supernova 
was a Typ e lib (Krause et al. I2OO8I: iRest et al.l 
200811201 ll) . 

Gonstraints on the global geometry of SN 
ejecta can also be obtained from polarization 
studies (see Wang fc Wheeler 20081 for a review) 
and from the line profiles in the nebular phase 



(e.g. Maeda et al. 


2008; 


Taubenberger et al.ll2009l: 


Milisavlievic et al 


2010 


). Both types of studies 



have revealed strong evidence for asymmetries, 
but the lack of spatial information means that 
there are no unique solutions for the 3-dimensional 
ejecta distributions. 

SN 1987A offers a unique opportunity to study 
the spatial distribution of ejecta from a core col- 
lapse SN. Its location in the LMG, only ^ 50 kpc 
away, together with the fact that the explosion 
was seen only 25 years ago (on 1987 February 23), 
make it possible to spatially resolve the ejecta and 
study the morphology in the homologous phase, 
before the inner ejecta crash into the circumstel- 
lar medium. 

The circumstellar material around SN 1987A is 
most apparent as three rings, inc lined by 38°-45 



with respect to the line of sight ([Tziamtzis et al 
201l[ ). The inner, equatorial ring spans ~ 
1.1" X 1.6" while the two outer rings, located 
above and below the equatorial plane, are about 
three times larger. The rings are thought to 
have been created some 20, 000 years before 
the exp losion, possibly as a result o f a binary 
merger (|Morris k Podsiadlowskill2007t) . By 1995, 
^ 2 900 days after the explosion, the outer- 
most part of the ejecta had reached the equa- 



1997 



called hot-spot to ap pear (jGarnavich et al 
Lawrence et al horn . The ring has since been ob- 



served to brighten exponential ly across the entire 



electromagnetic spectrum (e . g. iGroningsson et al 



2008a- Racusin et al.l 120091: IZanardo et al.l l2010t 



Dwek et al.ll2010l ). although there are now signs 

Dti- 

ll 



tha t the flux has started to level off in the opti- 
cal ( Groningsson et al.ll2008al : iMattila et al.. 201( 



and possibly also in soft X-rays (|Park et al.ll2011t 
Maggi et al.ll2012t ). 



To study the inner ejecta it is crucial to ob- 
tain high enough spatial resolution to separate it 
from the circumstellar rings. This can be achieved 
with the Hubble Space Telescope (HST) at opti- 
cal wavelengths and with the use of adaptive op- 
tics from the ground at near-infrared (IR) wave- 
lengths. We have recently used HST imaging 
observations, carried out approximately once per 
year since 1994 (day ^ 2, 800), to mea sure the op- 



tical light curve of the inner ejecta dLarsson et al 



201 ll ). This showed that the flux declined as 
expected from radioactive decay of ^^Ti until ~ 
5, 500 days after the explosion, after which the 
flux started to increase. The fact that '''^Ti is the 
dominant radioactive isotope at these late times 
is also supported by the observat ion of the asso- 



ciate d hard X-ray emission lines (|Grebenev et al 



20121 ). 

It is clear, however, that radioactive decay can- 
not explain the re-brightening of the ejecta. In- 
stead, we found that the most likely explanation 
for this is the energy input resulting from the X- 
ray emission produced by the collision between the 
outer ejecta and the equatorial ring. Because of 
the brightening it is possible to carry out a de- 
tailed study of the morphology of the ejecta even 
at very late times. 

There were early observations of asymme- 
tries in SN 1987A long before the spatially re- 
solved images of the ejecta were obtained with 
the HST. Between days 20 — 100, fine-structure 
appea red in the spectra (the so- called Bochum- 
event, iHanuschik fc Thimm I1990I and references 
therein), which was interpre ted as being due to a 



high-velocity clump of ^^Ni (lUtrobin et al. 



Early speckle obser vations (|Meikle et al 



19951) 



1987 



Nisenson et al. 



ments ( Jefferv 



198 



19911 



and polarimetry measure- 
also provided evidence for 
asymmetries. These early observations only probe 
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the structure of the outer hydrogen envelope and 
nearby gas though, which is not obviously con- 
nected to the inner ejecta. 

More recently, studies of the morp hology of the 



inner ejecta have been carried out bv lWang et al 



( |2002 ) , using HST sp ectral and imaging data from 
day ~ 4, 700, and bv lKiasr et all (|2010[) . using in- 
tegral field spectroscopy in the near-IR obtained 
with SINFONI at the V LT around day 6,800. 



While IWang et alj (|2002l ) favor a model in which 
the ejecta are prolate, with the long axi s pointing 
out of the equatorial plane, Kjasr et ahl (^2010) in- 
stead find that the ejecta are elongated closer to 
the plane of the ring. The latter conclusion is also 
supported by observations of light echoes, which 
probe the first few hundred days af ter the explo- 



sion from different viewing angles (|Sinnott et al 
20121 ). 



In this paper we study the morphology of the 
ejecta using HST imaging observations obtained 
between day 2,770 and 8,328, HST/STIS spec- 
tral observations from days 4,571, 6,355 and 8,378, 
as well as ground-based VLT/SINFONI integral 
field spectroscopy observations from days 6,816 
and 8,714. 

This paper is organized as follows: we describe 
the observations and data reduction in section 2, 
present an analysis of the images and spectra in 
section 3, and then use these results to explore 
some simple models of the ejecta in section 4. We 
finally discuss our results and present our conclu- 
sions in sections 5 and 6. 

2. Observations and data reduction 
2.1. HST Imaging 

HST imaging observations of SN 1987A have 
been carried out regularly since ^ 2, 700 days after 
the explosion (year 1994) as a part of the Super- 
nova 1987A INTensive Study (SAINTS). In order 
to study the temporal evolution of the morphology 
we selected six observations, separated by 1, 000 
days, taken in the F675W/F625W filters. These 
images are dominated by strong Ha-emission and 
offer the best signal-to-noise ratios. For two of the 
selected epochs (day 5,012 and day 8,328) there 
are also high-quality images in five other broad 
filters which we include in the analysis. The five 
filters together cover the entire wavelength inter- 
val between - 2, 000 and 10, 000 A, and thus 



allow a detailed study of the spectral energy dis- 
tribution. To show which lines contribute to the 
different filters we plot in Fig. [1] the filter through- 
put curves superposed on ejecta spectra from HST 
FOS-2 G570H (day 3,604, below ~ 6,800 A) and 
HST STIS G750L (day 4,381, above ~ 6, 800 A). 

We used the MULTIDRIZZLeEI software to re- 
move cosmic rays from the images, apply distor- 
tion corrections and combine in dividual exposures 
using the drizzle technique ([Fruchter fc Hook 
20021 ) . For each observation we drizzled the images 
to an optimal pixel scale based on the instrument 
and dither pattern used. The effective spatial res- 
olution of the resulting images is between 0.06" 
and 0.09". A set of common stars located around 
the remnant were used to align the images, giv- 
ing excellent accuracy with the positions of stars 
in the aligned images always agreeing to within 
0.02". The details of all the observations are sum- 
marized in Table [T] 

2.2. HST STIS observations 

The HST STIS observations that we use for this 
work were obtained on days 4,381 and 4,571 (year 
1999), day 6,355 (year 2004) and day 8,378 (year 
2010), using the G750M and G750L gratings. In 
the wavelength intervals covered by these grat- 
ings the Ha line dominates the emission from the 
ejecta. The spectral resolution is ^ 50 km s^^ for 
G750M and ~ 450 km s'^ for G750L. The details 
of the observations are summarized in Table [2] 

For each slit position there are at least three 
separate exposures taken at dithered positions 
along the slit. We combined these using the CAL- 
STIS software in order to remove cosmic rays and 
hot pixels. This successfully cleaned up the data 
from the first two observations, but a number of 
hot pixels remained in the observations from the 
later epochs. To remove these we used the IRAF 
fixpix task together with a user-supplied bad pixel 
mask. We then produced fiux-calibrated two- 
dimensional spectra from the cleaned data with 
the CALSTIS x2d task in the standard way, and 
finally extracted one-dimensional spectra by sum- 
ming over rows in the two-dimensional spectra. 
All spectra were corrected for t he 287 km s~^ sys- 



temat ic velocity of SN 1987 (jGroningsson et al 
2008al) 



^ http: / / stsdas.stsci.edu/multidrizzle 



3 



2.3. SINFONI observations 

We use integral field spectroscopy observa- 
tions of SN 1987A obtained with SINFONI at 
the VLT on days 6,816, 6824, 6825, 6839 and 
6843 (year 2005) and days 8,714 and 8,717 (year 
2011). For this work we consider only the obser- 
vations taken in the H-band, which is dominated 
by the 1.644 /im line. In the mod e hng o f the 



Kiasr et al.l (|2010t ) and 



"^'^Ti-powered sp e ctrurn , , ^ 

Ijerkstrand et al.l (|201l[ ) showed that this line is a 
blend of [Si I] and [Fe II], with [Si I] dominating, 
although both lines originate in the regions close 
to the inner Fe core. For the two observations that 
we consider in this paper the X-ray input may also 
add a contribution from the non-processed hydro- 
gen envelope. 



The SINFONI K band , discussed in lKiaer et al 



(|201fll ) las too low signal-to-noise to be useful for 
this analysis, but we note that the morphology 
seen in this band (dominated by the He I 2.058 /im 
line) is consistent with that seen in the H-band. 
The total H-band exposure times are 4,200 s and 
3,600 s for 2005 and 2011, respectively. 

A detailed analysis of the 2005 data ha s al- 
ready been presented by iKiaer et al. ()2010l) and 
we refer the reader to that paper for details re- 
garding the data reduction. The 2011 data were 
reduced in the same way, but with a more recent 
version of the pipeline (v2.2.3, dated 2010-08-12). 
Because there are no non-variable sources within 
the field of view of the SINFONI observations, we 
used VLT/NACO observations from day 7,201 and 
8,644 (Mattila et al., in prep.) to improve the 
flux calibration. We extrapolated the continuum 
level of the SINFONI spectra to account for the 
small differences in wavelength coverage between 
the two instruments, and further corrected for the 
differences in observing epoch by assuming that 
the flux from the ring evolves in the same way as 
the soft X-rays. The latter assumption is based 
on the expectation that the H-band flux is likely 
to correlate with the optical flux, which is ob- 
served to correlate with th e soft X-ray light curve 
(|Gr5ningsson et al.ll2008al ). Since there is no pub- 
lished light curve for the ring in the H-band there 
is a considerable uncertainty in this assumption. A 
conservative estimate is that the SINFONI fluxes 
derived in this way have uncertainties of ±20% 
for day 6,816 and ±15% for day 8,717. The pri- 



mary reason for the larger uncertainty in the for- 
mer case is the larger difference between the SIN- 
FONI and NACO observing epochs. As for STIS, 
all the spectra were corrected for the 287 km s"^ 
systematic velocity of SN 1987. 

3. Data analysis 
3.1. Imaging 

In Fig. [5] we show contour plots of the ejecta 
produced from the F675W/F625W images in the 
six selected epochs between day 2,770 and 8,328. 
In the left panel the contour levels are linearly 
spaced between the minimum and maximum sur- 
face brightness for each epoch, while in the right 
panel the same levels are used for all epochs. 
We note that these images are completely dom- 
inated by Ha-emission (see Fig. [Ij. The cen- 
tral part of the Ha line has a small contribution 
from Ca I A6572 (originating in the 0/Si and S/Si 
zones, Ijerkstrand et al. 2011) and the blue high- 
velocity wing has a small contribution from [O I] 
AA6300,6364. In addition, [Ca II] AA7291,7323 
contributes to the red wing of the F675W fllter, 
while Na I AA5890, 5896 influences the blue wing 
of the F625W filter. All these lines, however, only 
give a marginal contribution compared to Ha. 

The most important result from Fig. [2] is that 
there is a clear transition around 5,000 — 6,000 
days from a roughly elliptical shape to a more ir- 
regular morphology, where the brightest part of 
the ejecta is an arc-like structure that is offset 
from the center of the remnant. We note that the 
transition in morphology coincides in time with 
the re-brightening of the ejecta. The early images 
also exhibit an asymmetry in the north-east and 
south-west edges of the ejecta due to one of the 
outer rings of circumstella r material. The outer 
rings are fading with time (jTziamtzis et al.ll201ll ) 
and are not seen at all in the later images. 

A natural question to ask is whether the 
changes in morphology can be explained simply 
by the fact that the spatial resolution has in- 
creased with time, primarily due to the expansion 
of the ejecta but also because of the change of 
instruments and dithering strategies. To address 
this issue we first degraded the resolution of the 
later, higher-resolution images to match the first 
observation, which has the lowest resolution. The 
results are shown in the upper panel of Fig. (3] 
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Next we scaled the final image from day 8,328 
to allow a comparison with the previous images. 
As the ejecta are expanding homologously, this 
was achieved by increasing the pixel scale by a 
factor given by the ratio between 8,328 days and 
the number of days since the explosion for each 
of the other observations. The results are shown 
in the bottom panel of Fig. [3l It is clear that 
the change in morphology cannot be explained 
by the increased resolution offered by the expan- 
sion. It is thus likely that the morphology reflects 
the energy source powering the ejecta, i.e. that 
the early morphology (before ^ 5, 500 days) re- 
flects the radioactive input from '*'*Ti, while the 
late morphology is driven by the external X-ray 
illumination. 

In order to investigate the two morphological 
stages in more detail we selected one epoch from 
each stage (day 5,012 and 8,328, respectively) for 
which there is HST imaging observations in six 
broad filters covering the entire wavelength inter- 
val between - 2, 000 - 10, 000 A (Fig. [1]). These 
images are shown in Fig. |4l The images from day 
5,012 all display a similar elliptical morphology, 
while all the images from day 8,328 reveal a simi- 
lar edge-brightened morphology. The only notable 
exception is that the ejecta emission in the F225W 
and F841W images from day 8,328 appear slightly 
more extended than in the other images from the 
same epoch. 

To examine whether the spectral energy distri- 
bution has changed between day 5,012 and 8,328 
we also measured the flux of the ejecta in all of the 
images, the results of which are shown in the left 
panel of Fig. [5l The plotted fluxes are the total 
fluxes in each filter, i.e. the measured flux density 
multiplied by the rectangular bandwidth of the fil- 
ter and the ratio between the peak throughput and 
the throughput at the reference wavelength. We 
performed the measurements using elliptical aper- 
tures with a position angle of 15° and an axis ra- 
tio of 0.76, approximating the shape of the ejecta. 
The semi-major axis of the aperture was 0.25" for 
day 5,012 and 0.42" for day 8,328, chosen to ac- 
count for the expansion between the epochs and to 
include as much of the ejecta as possible, while still 
avoiding the equatorial ring. This semi-major axis 
corresponds to ^ 4, 300 km s^^ in the plane of the 
sky. The sky background was measured as an aver- 
age of five empty regions located around the rem- 



nant, and subtracted from the ejecta fluxes. We 
note that the sky background is low, contributing 
less than 5% of the flux in any of the measure- 
ments. 

A bigger source of uncertainty is the back- 
ground due to scattered light from the brightening 
equatorial ring. While this background is negligi- 
ble in the early epoch it is fairly high in the late 
one. To estimate the contribution from the ring 
at 8,328 days we constructed a model for it in the 
F625W band (where it is brightest) using the Tiny 
Tini HST poi nt-spread-funct i on m odeling tool, 
as explained in iLarsson et al.l ( 201 il l . Using the 
model we found that the scattered light from the 
ring contained within the ejecta aperture corre- 
sponds to 1.3 % of the flux measured in a 0.35" 
wide elliptical annulus centered on the ring itself. 
Subtraction of this scattered light reduces the to- 
tal flux in the ejecta aperture by about 27 % in 
this band. We corrected the ejecta fluxes in the 
other filters by measuring the ring flux in each fil- 
ter and assuming that the same fraction of this 
flux (1.3 %) always falls within the ejecta aper- 
ture. This correction reduced the ejecta fluxes by 
11 — 27 % for the di fferent filters. As discussed in 
Larsson et al.l (|201lh . we estimate that the calcu- 
lation of the amount of scattered light is accurate 
to within 15 % (i.e. the correction in the F625W 
band is 27 ± 4 % etc.). 

The right-hand panel of Fig. [S] shows the flux 
ratios between the two epochs, divided by the ra- 
tios expected due to the small differences between 
the filters used for the two epochs (cf . Fig. [Ij . The 
latter ratios were calculated using the SYNFHOtH 
package with the spectra in Fig. |l]as input. The 
plot shows that the flux has increased by factors 
of 2 — 3 in most fllters, with the exception of the 
F255W/F225W band where the increase is a fac- 
tor ^ 5. 

As a flnal point we compare the HST imaging 
results with images of the [Si I]/[Fe II] 1.644 ^m 
line from the SINFONI observations on day 6,816 
and 8,714. Contour plots of the SINFONI data 
covering the ±3, 000 km s~^ interval around the 
line are shown in Fig. |6l There are some clear 
differences compared to the HST F675W/F625W 
contour plots in Fig. [5J The peak brightness of 



^http: / /www. stsci.edu/hst/observatory/focus/TinyTim 

•^http: / /www. stsci.edu/instituto/softwaro_hardware/stsdas/synphot 
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the SINFONI images is closer to the center of the 
remnant, near the 'hole' seen in the HST images. 
The ejecta are also slightly more elongated in the 
north-south direction in the SINFONI images. 

These differences are even more apparent in 
Fig. [71 where we over-plot contours from the HST 
F625W image from day 8,328 on the SINFONI day 
8,714 contour plot. The images are well aligned 
(the equatorial ring has the same position and 
shape) and the difference in observing date of 
340 days only introduces a ~ 5 % difference in 
the size of the ejecta. The spatial resolution is 
slightly better in the HST image, but this is not 
enough to e.g. hide the hole in the SINFONI im- 
ages. In particular, if the hole was absent from 
the SINFONI images due to the broader wings of 
the PSFs, we would not expect the brightness to 
peak near the center of the hole. The SINFONI 
images are also consistent with the HST/WFC3 
F160W images and the H-band images obtained 
with VLT/NACO (Mattila et al., in prep.). 

We thus conclude that there are significant dif- 
ferences between the morphology of the 1.644 /im 
emission and the Ha emission, noting again that 
the latter is similar to the morphology seen in all 
broad-band filters in Fig. |4l 

3.2. Spectra 

In order to obtain information about the three- 
dimensional distribution of material we need to 
complement the imaging observations with spec- 
tra. Because of the homologous expansion of the 
ejecta, the measured Doppler shifts are directly 
proportional to the distance along the line of sight 
from the center of the explosion, making the spec- 
tral information very powerful. For this analysis 
we consider HST STIS spectra of the Ha line taken 
on days 4,571, 6,355 and 8,378, as well as SIN- 
FONI spectra of the [Si I]/[Fe II] 1.644 ^m line 
from days 6,816 and 8,714. We note that there is 
also a STIS observation from day 4,381 which was 
obtained with a wide, 0.5" slit covering all of the 
ejecta. The full ejecta spectrum from this obser- 
vation is shown in Fig. [1] but we do not consider 
it in this section since it is consistent with the 
spectrum from day 4,571, which has much better 
spatial and spectral resolution. 

The STIS observations from day 4,571 probe 
the radioactively dominated phase. The obser- 



vations were obtained with the G750M grating 
for three adjacent slit positions (slit width 0.1") 
tilted by about 30° from the north, thus covering 
most of the ejecta (see Fig. [5]). We extracted one- 
dimensional spectra in four 0.15" regions along 
each slit as marked in Fig. |8l The resulting spa- 
tially resolved profiles are shown in Fig. [9] and the 
total profile is shown in Fig. [101 Due to the excel- 
lent spectral resolution we can clearly separate the 
narrow lines from the outer ring ([N II] A6548, Ha 
and [N II] A6583) from the broad ejecta profile. 
The latter extends out to at least ±4000 km s~^. 
The spatially resolved profiles show a trend of 
more blueshifted emission in the north and more 
redshifted emission in the south. The emission is 
also slightly brighter in the north. This is reflected 
in the full profile in Fig. [TOl (black line), which is 
asymmetric with a somewhat stronger blue wing. 
This asymmetry may partly be due to the presence 
of dust (see further discussion in section [S.3.3p . 

Compared to this early observation the STIS 
observations from days 6,355 and 8,378 have 
poorer spectral and spatial resolution, as they 
were taken with the G750L grating and a slit width 
of 0.2". For day 6,355 three adjacent slit positions 
oriented in the north-south direction cover most 
of the ejecta, while on day 8,378 there is only one 
slit position covering the central part (see Fig. [Tl]) . 
As above, we extracted one-dimensional spectra 
in 0.15" steps along each of the slits (Fig. [Tl]) . 
The results for the two epochs are shown in Figs. 
[12] and [131 respectively, and the full spectrum 
from day 6,355 is shown in Fig. EQ\ The 2D STIS 
spectra are plotted in Fransson et al. (in prep.), 
showing both the Ha and the [Ca II] emission. 

The SINFONI observations were carried out 
461 and 629 days later than the two STIS observa- 
tions, which is close enough to make useful com- 
parisons (the offsets correspond to ^ 7% of the 
age of the supernova). We therefore include the 
SINFONI 1.644 ^m fine profiles from days 6,816 
and 8,714 in Fig[T2l and [131 respectively. In order 
to extract SINFONI spectra from the same regions 
as the STIS spectra we used the geometrical cen- 
ter of the circumstellar ring as a reference point, 
noting that the size and shape of the ring is very 
similar in the optical and near-IR. The extraction 
regions are shown superposed on the SINFONI 
images of the 1.644 /im line (±3000 km s~^) in 
Fig. 111! There are some uncertainties in the com- 
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parison with the STIS spectra owing to the dif- 
ferent point spread functions (the SINFONI PSF 
has an enhanced core, with broader faint wings), 
the positioning of the STIS sUts, the image qual- 
ity etc; however, these effects are smah compared 
to the systematic effects introduced by contami- 
nation from the circumstehar material. 

At these late epochs, the fading emission from 
the outer ring, which was identified in the STIS 
spectrum from day 4,571, will only make a small 
co ntribution to the e.iect a spectra (see light curves 
in lTziamtzis et al1l201l[ ). On the other hand, the 
equatorial ring is much brighter at these times 
(e.g.. Fig. m and it is likely to contribute to the 
spectra in Figs. [T^] and [T31 To get an estimate 
of the level of contamination we extracted spectra 
from the ring for both SINFONI and STIS and 
re-scaled them to the ejecta spectra, as shown in 
Fig. [H For SINFONI it is clear that the narrow 
(±250 km s^^) central component originates from 
the ring, but that there are no other strong lines 
from the ring that could contribute to the broad 
ejecta profile. In the case of the STIS observa- 
tions, it is harder to determine the contribution 
from the ring due to the lower spectral resolu- 
tion. The ring spectrum contains intermediate- 
veloci tv hues of [N 11] A6548, Ha and [N II] A6583 
(e.g. Groningsson et al. 2008b[ ). which together 
make a broad, skewed component that is likely 
to contaminate the ejecta spectra between about 
~ -1000 - -f 1800 km s"^ There are no strong 
emission lines from the ring that could contam- 
inate the broad wings of the ejecta Ha-profile 
though. 

Inspection of the spatially resolved ejecta pro- 
files in Figs. [12 and [T3] show the same trend of 
more blueshifted emission in the north and more 
redshifted emission in the south, as already noted 
for the earlier STIS observations in Fig. [S] The 
Ha and [Si I]/[Fe II] profiles mainly differ in their 
blueshifted emission, with the Ha-lines exhibit- 
ing smoother wings extending to higher velocities. 
These properties will be discussed in more detail 
in section [5] below. 

As the slit positions for the STIS observations 
on day 6,355 cover most of the ejecta we also 
show the full profile from these observations in 
Fig. [TOl together with the total profile from day 
4,571. Even though the spectra were extracted 
such that most of the ejecta would be included for 



both epochs, it is clear that some uncertainties are 
introduced by differences in gratings, slit width 
and orientation, as well as the expansion of the 
ejecta. Nevertheless, it is interesting to note that 
the two STIS profiles are very similar despite the 
change in morphology that has occurred during 
the five years between the two observations. The 
fluxes are also rather similar, which is likely due 
to a combination of the uncertainties mentioned 
above, together with the fact that the minimum 
flux was reached in the time interval between the 
observations 

The SINFONI data also cover all of the ejecta 
and we therefore next consider these observations 
in more detail. Compared to the STIS observa- 
tions, the SINFONI data have the advantage of 
higher spatial resolution perpendicular to the slit. 
In Fig. [15] we plot the three-dimensional emissivity 
in the [Si I]/[Fe II] line, obtained by converting the 
velocity along the line of sight to a distance from 
the center for each pixel, assuming homologous ex- 
pansion and a distance of 50 kpc to SN 1987A. The 
position of the equatorial ring and our line of sight 
are also included as reference points. Animations 
showing the f ull 3D view of the emissivity, create d 
with Mayavi ( Ramachandran fc Varoauauxll201l[) , 
are included in the online version of the journal. 

These plots and animations show that the emis- 
sion is concentrated in two "clumps", one to the 
north on the near side and one to the south on the 
far side with respect to the plane of the sky. In or- 
der to show the large-scale structure as clearly as 
possible we have chosen to only plot ejecta emis- 
sion brighter than 3 times the continuum level 
with velocities lower than 3, 500 km (excluding 
the central ±250 km s^^ along the line of sight). 
The first condition ensures that the plotted data 
points are well above the noise level and highlights 
the location of the brightest emission. We note 
that choosing a lower cutoff simply makes the ob- 
served structures slightly larger. The cuts in ve- 
locity were chosen in order to avoid emission from 
the ring blending with the ejecta emission. The 
downside of this is that some bright emission in the 
southern tip of the ejecta is excluded (cf. Fig. [6]) 
and that there is a small gap in the central region. 
We stress that none of these different cuts change 
the main conclusion that the emission is concen- 
trated to two clumps. It is also clear that there is 
very little change in morphology between the two 
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observations. 

From these three-dimensional maps we also cal- 
culate the radial intensity distribution by sum- 
ming over spherical shells of width 250 km s~^. 
The results are shown in Fig. [161 For both epochs 
the intensity peaks around 2, 300 km and the 
width of the distribution at half maximum is 
1, 700 km s~^. We note that there are systematic 
uncertainties in these distributions due to difficul- 
ties in defining the continuum level (< 15 % uncer- 
tainty) and the fact that the resolution is higher 
along the line of sight than perpendicular to it by 
a factor of 7 — 9 (the velocity resolution perpen- 
dicular to the line of sight improves with time due 
to the expansion). 

The change in the ejccta flux between the two 
SINFONI observations is a factor of 1.2 ± 0.3. 
While this is consistent with the HST observa- 
tions, the rather large uncertainties in the flux 
calibration (discussed in section 12. 3p means that 
we cannot exclude the possibility that a signifi- 
cant fraction of the emission at these wavelengths 
is still due to radioactivity. 

4. Modeling 

Because the ejecta look roughly elliptical on the 
sky in both the early HST images and in the SIN- 
FONI images it is interesting to investigate how 
well the spectra and images can be described by 
a model where the ejecta are distributed in the 
shape of an ellipsoid. We do not expect such a 
simple model to give a good fit to the data, but it 
is instructive to explore easily parameterized mod- 
els of ellipsoids with different shapes, orientations 
and dust opacities in order to gain insight about 
what kinds of geometries are possible or can be 
ru led out. The modeling done here improves that 
m iKiaer et all (|2010h in that we now include the 



effects of dust and also consider the spatially re- 
solved line profiles. 

In terms of comparing models with observa- 
tions we focus on those observations that have high 
enough spectral resolution to clearly distinguish 
the narrow lines from the equatorial and / or outer 
rings from the ejecta emission. This is possible for 
the STIS G750M observations from day 4,571 and 
for the SINFONI observations. As there is very 
little change between the two SINFONI observa- 
tions on days 6,816 and 8,714, we model only the 



first observation, where it is easier to separate the 
ejecta from the equatorial ring. 

In order to avoid degeneracies in the model we 
assume that the ejecta emission is uniform within 
the ellipsoid and that the distribution of dust is 
also uniform. We define the optical depth (r) as 
e""^ — (observed fiux from ellipsoid) / (total fiux 
emitted from ellipsoid). We further fix the posi- 
tion a ngle at the observed value of 15° (jKiaer et al 



20101) and require that the ellipsoid has the correct 



projected size on the sky for the two epochs. The 
latter condition means that we fix one axis of the 
ellipsoid (a) to the observed projected minor axis, 
and consider combinations of the other two axes 
{b and c > b) that give the correct observed pro- 
jected major axis. The parameters that we vary 
in the model are thus the inclination (9) of the 
longest axis, the combination of b and c, as well as 
the optical depth. 

Since our goal is simply to get an idea of possi- 
ble/impossible geometries we vary the parameters 
in rather big steps. We change the inclination in 
steps of 22.5°, consider eight allowed combinations 
of 6 and c for each inclination, and finally vary the 
optical depth between — 2.0 in steps of 0.5 for 
each set of 6, b and c. We evaluate the models by 
performing a fit to the ejecta profiles, excluding 
the narrow lines from the fit. 

None of the models manage to fully capture the 
observed profiles, thus confirming that the mor- 
phology is more complex than an ellipsoid. In 
Figs. [T8land[T9lwe show the expected line profiles 
for five examples of ellipsoids with different incli- 
nations and sizes, compared to the observational 
data from day 4,571 and day 6,816, respectively. 
All these examples have r = 0.5, since this gave 
the best overall match to the observed line profiles. 
This optical depth should only be seen as indica- 
tive, since it is clear that the ellipsoid models are 
too simple to explain the observations (see section 
15.3.31 for a discussion of dust in the ejecta) . 

It is clear from Figs. [T8l and fT9l that the models 
with the longest axis perpendicular to the equa- 
torial plane produce too much redshifted emission 
in the north and too much blueshifted emission in 
the south, while the models with the major axis 
parallel to the equatorial ring have the opposite 
problem. The problem is less severe if the axis ra- 
tios are lower, but in that case the model profiles 
are too wide in the central parts of the ejecta. 
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The model with the major axis in the plane 
of the sky actually describes the SINFONI spec- 
tra in the northern half of the ejecta reasonably 
well, but shows larger deviations from the obser- 
vations in the southern part. For the STIS data 
the same model has spectral peaks in approxi- 
mately the right places in the whole ejecta, but 
fails in that the wings of the central profiles are 
much sharper than the observations. We discuss 
the implication of these modeling results further 
in section [5] below. 

5. Discussion 

5.1. Thermalization of the radioactive and 
X-ray input 

Our observations provide us with intensity 
maps in different spectral lines. Because our main 
goal is to infer information about the spatial dis- 
tribution of the different elements, we need to 
discuss the connection between the energy input 
in the form of radioactivity and X-rays, the den- 
sity distribution, and the observed emis sion. This 
i s disc u ssed in detail in e^g^^ Kozma fc FranssonI 
( 19921 ). Ijerkstrand et all ( 2011 ) and Fransson et 
al. (2013, in prep. F13 from here on), so here we 
only summarize the main points. 

The de cay of '^'^Ti produces g amma-rays and 



positrons. iJerkstrand et al.l (j201ll ) estimate an op- 
tical depth to the gamma-rays of r-y w 0.05(^/8 yrs)" 
At 20 years is only 0.01. Wit h full trapping of 
the p ositrons, as indicated in Jerkstrand et al.l 
POIII ). almost all radioactive input will be de- 
posited in the Fe-rich material. The energy of the 
positrons is deposited as ionizations, excitations 
and heating, resulting in optical and IR emission 
lines. Because the energy input mainly occurs in 
the Fe-rich material, with a minor contribution 
from the Si/S zone, most of the resulting emission 
will come from ions of Fe-rich elements. 

In addition to the energy input directly con- 
nected to the ^''Ti decay, there is an important 
contribution from freeze-out emission in the hy- 
drogen envelope, where the low density makes 
recombination slow and energy input from ear- 
lier r adioactive decays of ^^C o and ^''Co domi- 
nate (jFransson fc Kozmalll993l ). This component 
mainly results in Lya, two-photon and Balmer line 
emission. 

Even in the hydrogen and helium zones the op- 



tical depth to the UV resonance lines from Mg I- 
II, Ti I-II, Fe I-II and other trace elements is large 
enough to scatter nearly all of the radiation emit- 
ted in the UV. Some of this will escape as blended 
lines at longer wavelengths in the UV and near- 
UV, while some is emitted as fluorescence radia- 
tion in the optical and IR. Spatially, this emission 
comes from the hydrogen envelope, even though 
the primary origin may be from the processed Fe 
core. 

At later epochs, when the ejecta are powered 
by the X-ray emission, the thermalization of the 
fast photoelectrons created by the X-ray absorp- 
tion is similar to that of the gamma-rays and 
positrons resulting from the radioactive decay. As 
long as the energy of the photoelectron is larger 
than ~ 0.3 keV, the fraction of the energy go- 
ing into heating, ionizations and excitations is in - 
dependent of energy (e.g., IXu fc McCravl Il99l[ ). 
It is, however, sensitive to the electron fraction, 
which determines the efficiency for Coulomb heat- 
ing. Estimates in F13 show that the ionization 
in the outer regions close to the ring can be high, 
> 0.1, while that in the core is typically ^ 10""^. 
For ionization fractions in this range, the total 
fraction of the X-ray energy emergin g as Ha is 
2-5% (|Xu fc McCravlll99ll : [Kozm a fc Fransson 



19921 ). Most of the remaining energy emerges as 
two-photon emission and Lya. In the same way as 
described above, this emission may be scattered by 
UV resonance lines. 

5.2. Evolution of the spectral energy dis- 
tribution 



In Larsson et al. 



(|201ll) we showed that the 
ejecta of SN 1987A have brightened by a factor 2-3 
in the F675W/F625W and F439W bands since day 
5,012. We also showed that the brightening can 
be explained by a change in the dominant energy 
source, from radioactive decay before day 5, 500 
(year 2001) to X-ray input from the ejecta - ring 
collision thereafter. We have now extended the 
analysis to a broader wavelength range, using five 
broad filters that cover the entire interval between 
~ 2, 000 and ~ 10, 000 A. As shown in Fig. U the 
flux has increased by a factor 2-3 between days 
5,012 and 8,323 in all the filters, with the excep- 
tion of the UV-band (F255W/F225W), where the 
increase is a factor of ^ 5. 

There are no spectra in the UV-band at late 
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epochs, so we cannot say exactly which hnes are 
responsible for the large increase, or if the flux 
correction between the F225W and F255W filters 
deduced from the early spectra has changed. How- 
ever, we note that at earlier epochs the band is 
dominated by Fe I-II and Mg I-II lines (Fig. [1]) 
and that most of the flux in these lines is due to 
resonance scattering in the hydr ogen-rich zones of 
the s upernova (section 15.11 and I Jerkstrand et al 
201 ih . As we will discuss below, most of the X- 



rays from the shocks are absorbed in the hydrogen- 
rich outer ejecta, so it is reasonable that the largest 
increase is seen in this wa velength interval. I n fact , 



from figures 3 and 4 in I Jerkstrand et al.l (|201lh 



it is clear that the core emits primarily in the 
3,000 - 6,000 A range, while the envelope dom- 
inates the emission at shorter and longer wave- 
lengths. This picture is also consistent with the 
fact that the filter with the second to largest 
increase is F814W, at the longest wavelengths 
(Fig.©. 

The evolution of the broad-band fluxes can 
be compared with the spectral evolution deter- 
mined from ground-based observations obtained 
with UVES at the VLT (F13). Although the spa- 
tial resolution is limited, the high spectral resolu- 
tion of the UVES data makes it possible to fol- 
low several lines from the ejecta. In particular, 
the broad emission lines from the inner ejecta can 
be separated from the narrow lines from the ring, 
as well as the boxy, broad lines from the reverse 
shock. Using this method, F13 found an increase 
in Ha by a factor of 4 — 6 between 5,000 and 8,000 
days, and an increase in the [Ca II] AA7291, 7323 
hnes (which fall within the F675W and F814W 
filters) b y a factor of 5 — 8 over the same period 



see also lLarsson et al.ll201lD . Both these changes 



are consistent with our HST broad-band photome- 
try within the observational and systematic errors, 
which include contributions from other lines, the 
continuum and scattered ring emission. A direct 
comparison of absolute flux levels determined from 
HST and UVES is not possible due to systematic 
differences in the contamination by emission from 
the ring, slit losses etc. 

The flux in the 1.644 fim line, which we mea- 
sured from the SINFONI observations on days 
6,816 and 8,328, has also increased, but the rather 
large uncertainties in the absolute flux calibration 
prevents us from making a meaningful comparison 



with the other filters. As already noted above, we 
cannot exclude the possibility that a significant 
fraction of the emission at these wavelengths is 
still due to radioactivity. 

5.3. Evolution of the morphology 

5.3.1. Edge-brightened morphology due to X-ray 
illumination? 

In this paper we have shown that the morphol- 
ogy of the ejecta seen in the HST images changes 
from an elliptical shape to a more irregular, edge- 
brightened geometry around the same time as the 
dominant energy source changes. A likely expla- 
nation for the late-time morphology is thus that 
it is due to the external X-ray illumination. The 
details of the X-ray deposition in the ejecta are 
discussed in some detail in F13, so here we only 
make some qualitative remarks that are relevant 
for the morphology. 

The location where X-rays from the shock re- 
gion are absorbed depend on the X-ray spectrum. 
Specifically, the column density penetrated by an 
X-ray photon with energy i? is oc iJ^/^. As a re- 
sult, the soft X-ray/extreme UV emission from the 
slow shocks in the ring will be absorbed in the out- 
ermost parts of the ejecta, while the hard emission 
(> 2 — 3 keV) can penetrate further in. 

As shown in IPark et aL ( 2011 ) there is a con- 
siderable and increasing hard flux from the shocks 
in the equatorial ring. Most of the hard emission 
is likely to be absorbed at the boundary between 
the envelope and core region, where the steep den- 
sity gradient and the dramatic increase in metal- 
l icity leads to an in crease in the X-ray opacity 
( Larsson et al. 201 1[ F13). A central low surface 
brightness region, largely protected from the X-ray 
illumination from the ring, is therefore a natural 
outcome, as also demonstrated by the modeling in 
F13. 

In view of this model it is also interesting to 
consider the top panel of Fig. [31 which shows the 
temporal evolution of the R-band morphology to- 
gether with a circle that follows the homologous 
expansion. In most of the ejecta the expansion 
follows the homologous flow. This is expected in 
both a radioactively- (or freeze-out-) dominated 
scenario, as well as in the case where most of the 
X-rays are thermalized in the steep density gradi- 
ent between the envelope and the core. However, 
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the southern part of the ejecta clearly 'expands' 
faster. This could indicate a high-density region 
at high velocities in this part of the ejecta, which 
was not reached by the early radioactive input but 
which is now being illuminated by the X-rays. Al- 
ternatively, this part of the ejecta is already being 
affected by the reverse shock. Both these scenarios 
point to an external energy source. 

Looking in more detail at the late-time mor- 
phology in the different bands in Fig. 21 we see 
that the emission appears somewhat more spa- 
tially extended in the F225W and F814W filters 
than in the other bands. In the former case this 
is likely due to the resonance scattering discussed 
in section 15.11 The forest of strong lines in the 
2000-3000 A region have large optical depths and 
may scatter emission at lower densities and thus 
larger radii than in the redder bands. 

In the case of the F814W band the emission 
is dominated by the Ca II triplet and the Mg 
II AA 9218.2,9244.3 lines (see discussion in F13). 
The Mg II emission is likely to be a result of Lya 
fluorescence in the hydrogen envelope, which may 
result in a more extended emissivity. However, 
confirmation of these issues will require a detailed 
spectral modeling, including the X-ray input. 

One possible explanation for the more centrally 
peaked morphology in the SINFONI images is 
that we are still probing the radioactive input 
with these observations. However, any hard X-ray 
emission penetrating the inner core of the ejecta 
would also result in stronger [Si I] and [Fe II] emis- 
sion from this region compared to Ha. The rea- 
son for this is that the non-thermal excitation by 
the electrons produced in the X-ray photoioniza- 
tion is roughly proportional to the abundance (in 
contrast to thermal excitation, which is mainly 
temperature-dependent) and the abundance of Si 
and Fe is higher in the core. 

In this context it is also important to point out 
why the morphology seen by SINFONI differs from 
the morphology in e.g. the F225W/F255W band, 
which is also dominated by metal lines. As dis- 
cussed above, the emission in this band is domi- 
nated by resonance scattering in the hydrogen-rich 
gas. A large fraction of the UV/optical emission 
from the metal core will thus have its last scatter- 
ing in the hydrogen envelope, while the emission 
at near-IR and longer wavelengths can escape di- 
rectly from the core. 



5.3.2. Correlation between ring and ejecta emis- 
sion 

If the ejecta are powered by radiation ema- 
nating from the circumstellar ring, we would ex- 
pect some degree of spatial correlation between 
the flux from the ring and the ejecta (we have 
already shown tha t the total fluxes correlate. 



Larsson et al.ll2Qlll ). The ~ 6 month light travel 
time between the ring and the inner ejecta is, 
however, too short to probe with the available 
observations. 

The first optical hot-spot ap peared i n the 
NE part of the r ing (jSonncb orn et al. IQQSt 



Lawrence et al.|[2000l ). and a corresponding bright- 
ening can indeed be seen in the closest part of the 
ejecta (days 5,012 and 6,122 in Fig. [2]). As more 
and more optical hot-spots appeared a larger frac- 
tion of the ejecta started brightening, as is ap- 
parent from Fig. [2] A caveat here is of course 
that the hot-spots observed in the optical/UV 
are produced b y radiative shocks with v elocities 



500 km s-i iGroningsson et al.l lioOSal ) . They 



< 

therefore mainly reflect the softest X-rays with 
energies < 0.3 keV, which are expected to be ab- 
sorbed in the outermost parts of the ejecta. The 
resulting effect on the morphology is difficult to 
judge due to the low surface brightness of this 
region. Slow recombination (freeze-out) may also 
delay the optical ejecta emission compared to the 
X-rays from the ring. 

For the harder X-rays, which are more im- 
portant for the deposition further in, we have 
to compare with the Chandra images. A set of 
0.3 - 8 keV images from day 4,71 1 to day 7,987 



are sh own in iPark et al.l (l2007|) and iRacusin et al 
()2009l) . Unfortunately, no separate images are 
shown for the different X-ray energies, although 
we know that the soft (0.5 - 2 keV) light curve 



(Park et al. 


201 if), but similarly to the opti 


Lg emission ( 


Groningsson et al. 2008a ) . Fur 



cal : 

thermore, Chandra observations at around 5,000 
days showed a different morphology above and be- 
low 1.2 keV, with the high-energy ban d exhibiting 



simil ar ities with the radio morphology ([Park et al 
2002t ). IPark et al.l (|2007t ). however, state that the 



later images up to 2007 show a similar morphology 
in the hard and soft X-ray bands. 

Although the resolution of the Chandra images 
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is considerably lower ('^ 0.5") than the HST im- 
ages, we note that the brightening in the X-rays 
first occurred in the NE and S E parts of the ring 
(see Fig. 1 of Park et al]|2007t ). In particular, the 
NE brightening corresponds to a similar brighten- 
ing of the ejecta in the HST image at 5,012 days 
(Fig. [2]). The SE X-ray brightening does not, 
however, have a prominent correspondence in the 
ejecta. We discuss possible reasons for this be- 
low. The next region to brighten in the Chandra 
images is the NW part, which increases dramati- 
cally between day 5,791 and day 6,533. Again, a 
similar evolution can be seen in the day 6,122 and 
day 7,226 ejecta images. The last Chandra region 
to turn on is the SW, which is also followed by 
a brightening in the HST ejecta images after day 
7,226. 

To summarize, we find a fairly good correspon- 
dence between the X-rays emission from the ring 
and the HST images of the ejecta, as expected 
from the X-ray dominated scenario for the en- 
ergy input. A more detailed analysis is beyond 
the scope of this paper, as it would have to take 
into account the difference in the azimuthal X-ray 
illumination between the soft and hard bands, the 
non-spherical distribution of the ejecta, as well as 
the possibility of dust obscuration. 

5.3.3. Effects of the Ni-bubble and dust 

Although it is clear that the X-ray illumination 
can explain the edge-brightened morphology in the 
HST images, it is also possible that the observed 
morphology is affected by the 'Ni-bubble' or the 
presence of dust in the ejecta. We now discuss 
these alternatives in turn. 

The Ni-bubble is expected to occur as a re- 
sult of heating by the ^^Ni de cay during the 
first days after the explosion (e.g.,|l 



1992t iLi et al.l [l99l lEUinger et al 



erant fc Benz 



20121 ). The 



bubble ( which may be composed of many small 
bubbles, iLi et al.l 19931) has the effect of pushing 
the nuclear burning zones located outside the iron 
core outwards, thus creating a dense shell sepa- 
rating the hydrogen- to oxygen-rich zones from 
the iron-rich material. This could in principle ex- 
plain the difference in morphology between the 
SINFONl and HST images (Fig. ^. 

However, we believe that this explanation 
is unlikely since we find that the ejecta mor- 



phology has changed considerably over the ob- 
served time interval, in contradiction with the 
expectation that the Ni-bubble should be frozen 
into the homologous expansion. An additional 
problem with this explanation is that the fill- 
ing factor of t he Ni-bubble was estini ated to 
be only ~ 0.2 ( Kozma fc FranssonI 1998 ) within 
~ 1,800 km s~^, which is likely to be too small 
to explain t he observations. Th e fillin g factor es- 
timated by iKozma fc FranssonI (Il998l) is smaller 



Li et al.l (|l993l) {who find 



than the results of 
filling factor of ~ 0.6), but the latter work as- 
sumed that no primordial Fe contributes to the 
line emission. 

There is plenty of evidence for the formation 
of dust in SN 1 987A, both from line profiles 
(jLucv et al ] Il989|)_ and thermal emission during 



the first years (IWooden et al.l Il993f). as well as 



recent Herschel ( Matsuura et al .1 l201l[ ) and sub 



millimeter ( Lakicevic et al. I I2OIII) observations, al- 
though the amount of dust formed is highly de- 
bated. Any dust present in the ejecta is likely to 
have formed in the metal-rich core of the SN, and 
could thus cause an obscuration of the central re- 
gions. The difference between the SINFONl and 
HST images could then partly be an effect of the 
lower extinction in the near-IR SINFONl band. 

In order for this effect to be important the dust 
has to be optically thick in the F814W filter at 
~ 0.85 /^m (where the hole is seen. Fig. S]), but not 
in the SINFONl H-band at ~ 1.6 /xm. The com- 
paratively small difference in extinction between 
0.85 /im and 1.6 /im (a factor of ^ 3) then re- 
quires that we observe at just the right epoch when 
these conditions are met. Since the optical depth 
of the dust is expected to decrease as where t 
is the time since explosion, substantial fine tuning 
is required in order for optical depth unity to be 
reached at ~ 20 years. In addition, the decrease 
in optical depth would be considerable from 2005 
to 2010, leading to an increased transparency at 
Ha, which is not indicated in the HST images. 

It is more likely that the dust is in the form 
of very optically thick clumps, which would give 
the same extinction at optical and infrared wave- 
lengths. This is supported by the fact that the [O 
I] AA 6300, 6364 lines, which became asymmetric 
after about 500 days, presumabl y as a result of 
dust formation ( Lucv et "al1ll989[ ) still have very 
similar profiles around 4,000 days (F13). These 
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l ine profiles, as well a s the Mg I] A 4571 line 



(jJerkstrand et al.l 120111 ) indicate an effective op- 
tical depth close to unity, which in the case of op- 
tically thick clumps represents a covering fraction 
of the clumps. 

Absorption by dust in the ejecta will primarily 
have the effect of reducing the redshifted emission. 
As clearly seen in Fig. [121 the red halves of the pro- 
files decay smoothly (and in a similar way for Ha 
and [Si I]/[Fe II]) while the blue profiles rise more 
sharply for [Si I]/[Fe II]. This indicates that Si/Fe 
is distributed in more sharply defined clumps than 
H, but that both components are affected by dust 
absorption. More definite conclusions regarding 
the effects of dust on the observed morphology 
and line profiles will have to await spatially re- 
solved images of the dust emission in the ejecta. 
This will hopefully be achieved with ALMA in the 
near future. 

5.4. Evolution of the line profiles 

Our study of the ejecta spectra at different 
epochs does not reveal any strong temporal evo- 
lution of the line profiles. In particular, we do 
not see any major changes between the Ha pro- 
files from day 4,571 (when radioactive decay was 
still dominating the energy input. Figs. l9l and [T0|) 
and day 6,355 (when the X-ray input had started 
to dominate, Figs. [T^] and [TUl). We note, how- 
ever, that the comparison is made difficult by the 
change of grating, slit width and slit orientation. 

The change in the dominant energy source is 
not expected to change the spectra dramatically. 
Both the radioactive decay and the X-ray emission 
result in secondary non-thermal electrons, which 
subsequently give rise to line emission via the same 
physical processes, thus producing similar sp ec- 
tra. (see section lSTI above. Larsson et al.l[2011 and 
F13.). In addition, the freeze-out effect in the en- 
velope will delay changes in the energy input. The 
main difference between the two phases is that the 
energy deposition occurs at higher ejecta velocities 
in the X-ray dominated phase. It is therefore rea- 
sonable to expect some changes to the line profiles 
at later epochs, since the location of the X-ray de- 
position is likely to change as the ejecta expand 
and the X-ray light curve evolves. 

No major changes have been seen yet, as evi- 
dent from a comparison of the spatially resolved 



profiles in Figs. [12] and [131 which show the Ha 
and [Si I]/[Fe II] emission around 6, 600 and 8, 500 
days, respectively. However, the one slit that was 
used in the latter STIS observation does not con- 
tain the blob in the western part of the ejecta, 
where the strongest brightening is seen in the HST 
images. 

An important lesson from this work is that the 
total line profiles of the ejecta are only slightly 
asymmetric (Fig. [Tn]), despite the large asymme- 
tries seen in the images and in the spatially re- 
solved line profiles. The fact that such prominent 
asymmetries can nearly cancel out in the inte- 
grated profile highlights the uncertainties involved 
in trying to obtain information about the explo- 
sion geometry from the spectra of distant super- 
novae, where no spatial information is available. 

5.5. Spatial distribution of the ejecta and 
implications for explosion models 

All the observations that we have considered 
in this paper show that the ejecta morphology 
is highly non-spherical. There is no well-defined 
symmetry axis for all the emission, but it is clear 
from Fig. [15] that the [Si I]/[Fe II] emission is closer 
to the plane of the ri ng than perpend icular to it 
(as already noted by Kiaer et al]l201Cl[ ). The Ha 
emission is consistent with having the same distri- 
bution relative to the ring plane, but the emission 
extends to higher velocities than [Si I]/[Fe II]. The 
fact that the ejecta are preferentially distributed 
close to the plane of the ring has also been inferred 



from observations of light echoes (jSinnott et al 



2012| ). although it should be noted that these ob- 
servations primarily probe the outermost ejecta. 

We also find that there is a north-south asym- 
metry for both [Si I]/[Fe II] and Ha; the strongest 
emission in the southern part is concentrated 
closer to the plane of the sky, while the strongest 
emission in the northern part is closer to the plane 
of the ring. As discussed above, this asymmetry 
may be partially enhanced by dust absorption. 

One of the main motivations for this study was 
to obtain information about the explosion geome- 
try. As was discussed in section [STTl the observed 
emissivity distribution does not necessarily trace 
the matter distribution, but is a complex convolu- 
tion of the input from radioactivity or X-rays, the 
density and abundance distribution, dust absorp- 
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tion, as well as the thermalization process. 

In the case of [Si I]/[Fe II] we are clearly seeing 
the innermost regions of the ejecta, and it there- 
fore seems hkely that the observed emission is a 
reasonably good representation of the density dis- 
tribution, regardless of what the dominant energy 
source is (see section [575)1 . In the case of hydro- 
gen it is harder to relate the observed morphology 
to the density distribution since the morphology 
changes with time, responding to changes in the 
energy source. However, for the X-ray dominated 
scenario we have argued that the hole and the ring- 
like structure can be interpreted in terms of the 
high-density core and the hydrogen envelope just 
outside the core, respectively. This means that 
the late-time observations still provide useful in- 
formation about the structure of the ejecta, even 
though the emission is not directly proportional to 
the density. 

Dust absorption is unlikely to significantly af- 
fect the observed morphology in other ways than 
reducing the intensity of the redshifted emission. 
In particular, we have shown in section U that an 
ellipsoid elongated perpendicular to the equatorial 
ring (i.e. the simplest form of a jet model) is in- 
consistent with the observations even when dust 
is included in the model. 

Keeping the above caveats in mind we com- 
pare the observed morphology with existing mul- 
tidimensional simulations of the explosion geome- 
try. In the absence of detailed predictions for the 
Ha and [Si I]/[Fe II] intensity distributions from 
these models we limit ourselves to some qualita- 
tive comparisons. The most detailed calculations 
of the explosion geomet ry are the 3D simulations 
bv lHammer et al. ( 2010l) . These calculations start 
from a self-consistent explosion model which pro- 
vides the initial non-spherical irregularities that 
act as seeds for the continued evolution. Unfortu- 
nately, the simulations do not include the heating 
by the ^^Ni decay, which creates the 'Ni-bubble'. 
They are also based on a 15 Mq progenitor, which 
is at the lower e nd of the published mass es timates 
for SN 1987A. iFransson fc Kozmal (|2002[ ) find a 
progenitor mass of 18 — 2 Mtn based on the esti- 
mated oxygen mass, while ISmartt (l2009l ) estimate 
a progenitor mass in the range 14 — 20 Mq- 

Compared t o simil ar, previous calculations in 



ity mixing of the different elements, including Ni- 
blobs with high velocities. These fast clumps 
have also been seen in later SPH simulations by 
E Uinger et al.l (|2012l) . Even though lEllinger et al 



torn do not start from a self-consistent explo- 
sion model, they do include the ^^Ni decay, and 
also run their simulations for a much longer time 
than Hammer et al. ( 2010l ). making sure that the 
effects of the radioactive heating are included. 

Starting with the [Si I]/[Fe II] emission, we note 
two prominent blobs in the 3D maps shown in Fig. 
[TSl The reduction in flux at low velocities seen 
in this plot is somewhat exaggerated because the 
±250 km s~^ interval along the line of sight was 
removed due to the strong contamination from the 
ring. However, a real decrease at low velocities is 
present, clearly seen in e.g. the dips on both sides 
of the narrow line from the ring in the northern 
part of the ejecta (Figs. [T2landfT3l). These c lumps 



can be compared with the iHammer et al. I (l2010h 
simulations (their Figs. 2 and 3), where some 
prominent 'Ni' clumps (which include both Si and 
Fe) extend out to high velocities. However, de- 
spite the qualitative agreement, our observations 
show a significantly higher degree of asymmetry 
on large scales. Our integrated radial distribu- 
tion of emission (Fig. [T6)) also peaks at a higher 
velocity (^ 2,300 km s~^) compar ed to the radial 



mass distribution of 'Ni' in Fig. 6 of lHammer et al 
(l2010[ )r~ 1,200 km s'^). 



The observations of Ha in SN 1987A have lower 
spectral resolution and are also harder to relate to 
the full density distribution. We therefore limit 
our comparison with explosion models to two very 
general conclusions that can be drawn from the 
observations. These are that the H emission ex- 
tends to higher velocity than Si/Fe and that the 
the emission decays more smoothly with velocity 
than Si/Fe. Both of these conclusi ons are in quali- 
tative agreement with the results of lHammer et all 
()2010l) . 



Kifonidis et al.l (|2006l ). these 3D calculations are 



found to give considerably more efficient veloc- 



From the late-time HST images we can also es- 
timate the size scale of different structures in the 
ejecta, bearing in mind that this only gives pro- 
jected sizes and that we cannot resolve structures 
smaller than ^ 500 km s^^. The diameter of the 
'hole' is around 1,000 - 2,000 km s'^ for all the 
late-time images, but the center of the hole is offset 
by ~ 800 km s~^ from the center of the remnant. 
The width of the ring-like structure corresponds to 
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a velocity width between ^ 1, 000 — 3, 500 km s~^. 
As noted above, we may identify the hole with the 
high-density core and the ring structure with the 
hydrogen envelope just outside the core. The peak 
of the emission in the SINFONI images is close to 
the position of the hole in the HST images, in 
agreement with this interpretation. 

In terms of studies of other young supernova 
remnants the most detailed results regarding the 
ejecta str ucture have been obtained in the case 
of Cas A (iFbs en et al.| j2006t iDeLanev et al.|[2010t 
Hwang fc Lamina, 20l3 ). A comparison with our 
results is not straightforward since both the su- 
pernova type and the type of observational data 
available are different. In par ticular, the fact that 



Cas A was a type lib SN (|Krause et al.l 120081 ) 



means that the H mass was very small; some- 
thing which will strongly affect both expansio n 



velocities and instabilities ( Iwamoto et al.l 119971 ) 



Compared to SN 1987A, the Cas A observations 
have dramatically higher spatial resolution, and 
also differ from our observations in that most of 
the observed emission comes from ejecta passing 
the reverse shock. Despite these differences, it is 
interesting to note that the knots, jets and pistons 
in Cas A ha ve been su ggested to have a flattened 
distribution (jPeLanev e t al. 2010,) . While we do 
not detect the same type of structures, we find 
that the ejecta in SN 1987A also have a fiattened 
distribution. 

In summary, we find reasonable qualitative 
agreement with the models, but we stress that 
a more detailed comparison requires both more 
realistic spectral modeling, as well as explosion 
simulations with parameters matching SN 1987A 
and all physical effects included. In terms of obser- 
vations, future ALMA observations of the spatial 
distribution of dust and molecules will help reduce 
uncertainties about the three-dimensional ejecta 
distribution. HST STIS observations with narrow 
slits covering the full ejecta would also greatly 
improve the comparison between the [Si I]/[Fe II] 
and Ha emission. 

6. Conclusions 

We have studied the morphology of the ejecta 
of SN 1987A using HST images obtained between 
day 2,770 and 8,328. In addition, we have ana- 
lyzed HST STIS spectra from three epochs after 



4,500 days, as well as SINFONI integral field spec- 
troscopy for two epochs after 6,800 days. The HST 
imaging covers the 2, 000—10, 000 A interval, while 
the STIS spectra mainly probe the Ha emission. 
For the SINFONI data we focused on the 1.644 ^m 
[Si I]/[Fe II] emission. Our main conclusions are 
as follows: 

• There is a strong temporal evolution of the 
morphology in the HST images, from a 
roughly elliptical shape before ~ 5, 000 days, 
to a more irregular, edge-brightened mor- 
phology (with a 'hole' in the middle) there- 
after. The transition is most likely due to 
the fact that X-ray emission from the equa- 
torial ring starts to dominate the energy in- 
put to the ejecta at this point. We find that 
other effects, from either the Ni-bubble or 
dust, are unlikely to significantly effect the 
morphology. The X-ray input also results 
in a brightening of t he ejecta, as we hav e 
previously shown in iLarsson et al.l (|201ll ). 
We note that the change in morphology and 
the late-time brightening is observed in all 
five broad HST filters between 2,000 and 
10,000 A, with the largest increase in flux 
occurring at the shortest wavelengths. 

• The SINFONI images of the 1.644 ^m [Si 
I]/[Fe II] emission show a more uniform mor- 
phology, with the brightest regions appear- 
ing close to the faintest, central regions of 
the HST images. This may be due to a 
remaining significant contribution from ra- 
dioactive decay in the inner part of the 
ejecta. However, because of the high abun- 
dances of Si and Fe in the core, a similar 
morphology would also be expected if there 
is a significant hard X-ray fiux reaching the 
innermost ejecta. Accurate flux measure- 
ments at these wavelengths over the coming 
years are needed in order to determine what 
the dominant energy source is. 

• The STIS Ha and SINFONI 1.644 ^m [Si 
I] / [Fe II] line profiles arc similar in the sense 
that they show predominantly blueshifted 
emission in the north and predominantly 
redshifted emission in the south, showing 
that the ejecta have a fiattened distribution 
fairly close to the plane of the equatorial 
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ring. There is also a north-south asymme- 
try: the emission in the northern half of the 
ejecta is concentrated around the plane of 
the ring, while the emission in the south- 
ern part is mostly above the ring, closer to 
the plane of the sky. This asymmetry may 
partly be due to dust absorption in the far 
side of the ejecta though, as suggested by 
the smooth decrease of the redshifted emis- 
sion. The Ha and 1.644 /xm [Si I]/[Fe II] 
profiles differ in that the Ha emission ex- 
tends to higher velocities, as expected from 
the differences between the images. 

We were unable to reproduce the observed 
images and spectra using a model for the 
ejecta in the form of a uniform ellipsoid 
with a uniform distribution of dust. The 
fact that the real distribution of ejecta is 
much more complex is also evident from the 
SINFONI 3D emissivity maps. Most of the 
SINFONI [Si I]/[Fe II] emission is concen- 
trated to two clumps in the velocity inter- 
val 1, 500 - 3, 000 km s'^. We do not have 
the spatial or spectral resolution to esti- 
mate the size of clumps from the Ha spec- 
tra, but in the HST images there is sub- 
structure with projected velocity widths of 
^ 800 km s"^ (the diameter of the 'hole') 
to ~ 3, 500 km s~^ (the thickest part of the 
ring structure). 

We have compared our results with the ex - 
plosion models by iHammer et al. (I2OIOI) . 
While there is some qualitative agreement, 
we also found that the observations show a 
higher degree of large-scale asymmetry than 
the simulation. We stress, however, that 
this explosion model was not set up specifi- 
cally to model SN 1987A and that radioac- 
tive heating was not included. In terms of 
explosion models we also note that the ob- 
served distribution of emission is inconsis- 
tent with a jet-like structure perpendicular 
to the equatorial ring. However, this does 
not necessarily rule out explosion models 
where most of the energy (but not ejecta) 
comes out along the poles. 
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Fig. 1. — System throughput for the HST filters used in our analysis, normalized to the peak value for each 
instrument. WFPC2 and WFC3 filters are shown as solid and dashed lines, respectively. Spectra of the 
ejecta from day 3,604 (HST/FOS-2/G570H, light gray) and 4,381 (HST/STIS/G750L, dark gray) are also 
shown, with the most important lines identified. 
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Fig. 2. — Contour plots of the ejecta produced from HST F675W/F625W images. The black star marks 
the position of the peak brightness in the first observation on day 2,770. The contours in the time-series 
to the left are linearly spaced between the minimum (blue) and maximum (red) for each epoch, while the 
images to the right have identical contour levels for all epochs. The red contours on white background 
correspond to higher surface brightness than the filled contours (see intensity scale to the right; the numbers 
are proportional to the flux per pixel). To demonstrate that the ejecta are spatially resolved in the first 
observation we show a nearby star, scaled to the same peak brightness as the ejecta, in the inset in the upper, 
left panel. The field of view for each of the images is 0.9" x 0.9". The first hot spot on the ring appeared in 
the north-east (outside the region shown here), i.e. in the same direction as the ejecta start re-brightening. 
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Fig. 3. — The top row shows contour plots of the ejecta as in Fig. [51 but all scaled to the spatial resolution 
of the lowest resolution images. The bottom panel shows the low-resolution version of the image from day 
8,328, reduced in size and "re-observed" at each epoch, assuming that the ejecta were expanding freely. The 
contour levels are linearly spaced between the minimum (blue) and maximum (red) for each image in order 
to highlight changes in morphology, and the white circle is included as a reference point that follows the 
homologous expansion (with 2, 800 km s~^ in the plane of the sky). The comparison clearly reveals that the 
change in morphology cannot be explained as a result of increased resolution due to the expansion. 




Fig. 4. — HST images from day 5,012 (left panel, taken with WFPC2) and day 8,328 (right panel, taken with 
WFC3) in six different optical filters. The field of view of each image is 1.9" x 1.5". The filter throughput 
curves are shown in Fig. [TJ Note that there are some differences in wavelength coverage for especially the 
F255W/F225W and F675W/F625W filters. The minimum and maximum levels in the images were set to 
highlight changes in the morphology. While they change from filter to filter, the scales are identical per filter 
for the two epochs. The bright point source in the bottom, right corner of the ring is a star. The bright 
point source in the upper left is the first hot spot appearing on the ring. 
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Fig. 5. — The left panel shows ejecta fluxes on day 5,012 (black circles) and day 8,328 (red triangles) for all 
the images plotted in Fig.lH The fluxes from day 8,328 have been corrected for scattered light from the ring. 
The errors are about the same size as the plot symbols . The fluxes have b een de-reddened by Eb_v = 0.19 
mag using the Cardelli et al. ( 1989( ) reddening law (see France et al. ( 2011 ) for a motivation of this choice). 
The right panel shows the ratio of the fluxes, corrected for the fact that the filters used for the two epochs 
have slightly different wavelength coverage and throughputs (cf. Fig.[T]). 
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Fig. 6. — Contour plots of the ejecta produced from SINFONI images of the [Si I]/[Fe II] 1.644 /im line 
in the ±3, 000km s~^ range. The contours are linearly spaced between the minimum (blue) and maximum 
(red) for each epoch. The black star marks the position of the peak brightness of the first HST image from 
day 2,770 (cf. Fig. [2]). The field of view is 0.9" x 0.9". 
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Fig. 7. — Comparison of the SINFONI [Si I]/[Fe II] image from day 8,714 (filled colored contours) and the 
HST F625W image from day 8,328 (black contours). The SINFONI contours are the same as in the right 
panel of Fig. [SJ while the HST contours correspond to every other level above the mid-point in the last image 
of the left time-series in Fig. [2l The field of view is 0.9" x 0.9". The peak brightness of the two images are 
clearly offset, with the SINFONI emission peaking closer to the center of the remnant, in the 'hole' of the 
HST image. Note also that the brightest region in the southern part of the SINFONI image coincides with 
the minimum in the southern part of the HST image (see also Figs. [2] and S]). 
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Fig. 8. — The three slit positions used for the STIS observations on day 4,571 are shown as dashed hnes 
superposed on the WFPC2 F675W image from day 4,727. The sht width is 0.1" and the field of view is 
2.2" X 1.4". The extraction regions used for producing the spectra in Fig. [9] are shown as solid rectangles 
(size 0.1" X 0.15"). The 0.5" slit that was used for the STIS observations on day 4,381 (spectrum shown in 
Fig. [T]) has nearly the same orientation as the 0.1" slits plotted here. 
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Fig. 9. — Spatially resolved STIS G750M spectra from day 4,571. The spectra are centered on zero velocity 
for Ha. Each spectrum was extracted from a 0.1" x 0.15" region, as shown in Fig. [51 The narrow lines 
originate from one of the outer rings, which is partly projected on top of the ejecta. 
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Fig. 10. — Comparison of ejecta Ha profiles froni HST STIS observations obtained on day 4,571 (G750M 
grating, black line) and day 6,355 (G750L grating, red line). The spectra were extracted from the total 
regions marked with solid lines in Fig. [5] and the top, left panel of Fig. [TI] respectively. The extraction 
regions cover most of the bright, central ejecta for both epochs, although the differences in resolution (nearly 
a factor 10 better in the early observation), slit width and orientation, as well as the expansion of the ejecta, 
introduce some uncertainties in the comparison of the full profiles. 
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Fig. 11. — Extraction regions (solid red lines) used for producing the spectra in Figs. [12] (top row) and fT3l 
(bottom row). All the extraction regions are 0.2" x 0.15" and the field of view for each of the images is 
1.8" X 1.1". Top left: The three slit positions used for the STIS observations on day 6,355 are shown as 
dashed lines superposed on the ACS F625W image from day 6,790. Top right: SINFONI image from day 
6,816 in the wavelength range corresponding to 1.644 /im ±3000 km s^^. Bottom left: The slit position used 
for the STIS observations on day 8,378 is shown as dashed lines superposed on the WFC3 F625W image 
from day 8,328. Bottom right: SINFONI image from day 8,714 in the wavelength range corresponding to 
1.644 ^m ±3000 km s^^. 
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Fig. 12.— Spatially resolved STIS G750L (red) and SINFONI H-band (black) spectra from day 6,355 and 
6,816, respectively. The STIS spectra are centered on zero velocity for Ha while the SINFONI spectra are 
centered on zero velocity for the [Si I]/[Fe II] line at 1.644 fim. The SINFONI spectra have been rescaled by 
a constant to match the flux level of the STIS spectra. Each spectrum was extracted from a 0.2" x 0.15" 
region, as shown in Fig. 1111 
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Fig. 13.— Spatially resolved STIS G750L (red) and SINFONI H-band (black) spectra from day 8,378 and 
8,714, respectively. The STIS spectra are centered on zero velocity for Ha while the SINFONI spectra are 
centered on zero velocity for the [Si I]/[Fe II] line at 1.644 fim. The SINFONI spectra have been rescaled by 
a constant to match the flux level of the STIS spectra. Each spectrum was extracted from a 0.2" x 0.15" 
region, as shown in Fig.HH 
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Fig. 14. — Ejecta spectra from the top, middle panel of Fig. [12] shown together with rescaled spectra of 
the circumstellar ring. STIS spectra (day 6,355, centered on zero velocity for Ha) are shown to the left 
and SINFONI spectra (day 6,816, centered on zero velocity for the 1.644 fim [Si I]/[Fe II] line) to the right. 
In both plots the scaled ring emission is shown as a dashed blue line, illustrating the maximum possible 
contamination from the ring. 
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Fig. 15. — Three-dimensional view of the 1.644 emission from the ejecta observed by SINFONI on day 
6,816 (top row) and 8,714 (bottom row) for different viewing angles. The blue sphere and line show the 
position of the observer and the line of sight. The black circle shows the peak of the emission from the 
equatorial ring, and the black line is perpendicular to the plane of the ring. The radius of the ring is 0.82", 
which corresponds to 6.1 x 10^*" cm at a distance of 50 kpc. The plots show ejecta emission brighter than 3 
times the continuum level for each epoch, with darker and bigger spheres representing brighter emission. The 
ejecta emission is plotted out to ±3,500 km s~^ (assuming homologous expansion), excluding the interval 
between ±250 km s~^ along the line of sight, which is contaminated by emission from the ring. Animations 
of the 3D emissivity is available in the online journal. 
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Fig. 16. — Radial distribution of the [Si I]/[Fe II] 1.644 /um emission from the SINFONI observations on day 
6,816 (red dashed hne) and 8,714 (black line), summed over spherical bins with velocity width 250 km s~^. 
The continuum was calculated over the 3, 500 — 4, 000 km s~^ interval and subtracted from the other bins. 
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Fig. 17. — Ellipsoid geometries used for producing the model spectra shown in Figs. [TBI and [T9l The size 
of the ellipsoids shown here give the correct projected size on the sky for the SINFONI observations on day 
6,816. The axis ratios of the red and blue ellipses are 2.2:1:4.2 and 1:1.2:1.5 while the cyan one has 1.2:1:1.6, 
where a:b:c are defined in the text. 
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Fig. 18. — Model Ha profiles for different ejecta geometries shown together with the spatiahy resolved 
spectra from the STIS G750M observations obtained on day 4,571 (Fig. [9]), extracted from the regions 
shown in Fig. \8\ The observed spectra have been binned by a factor of five compared to Fig. |9] and the 
narrow lines from the outer ring have been removed. All the ejecta models assume a uniformly emitting 
ellipsoid with a uniform distribution of dust and optical depth t w 0.5. For the red, blue and cyan lines the 
major axis of the ellipsoid is perpendicular to the ring plane, parallel to the ring plane and in the plane of 
the sky, respectively (cf . Fig. [T7| . The solid and dashed lines correspond to the more and less elongated 
ellipsoids in Fig. [T71 respectively. All these model geometries give the right observed shape as projected on 
the sky, but it is clear that none of them successfully describe all the spectral profiles. 
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Fig. 19. — Model [Si I]/[Fe II] profiles for different ejecta geometries shown together with the spatially 
resolved spectra from the SINFONI observations obtained on day 6,816 (Fig. [T2|) . extracted from the regions 
shown in Fig. 1111 The narrow spectral component from the equatorial ring has been removed from the 
spectra. All the model profiles are as described in the caption of Fig. [181 Note that the cyan lines provide 
an acceptable description of the spectral profiles in the northern half of the ejecta for these observations. 
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Table 1 
HST Imaging observations 



Date Days since explosion^ Instrument Filter Exposure time 



1994-09-24 


2,770 


WFPC2 


F675W 


600 


1998-02-06 


4,001 


WFPC2 


F675W 


400 


2000-11-13 


5,012 


WFPC2 


F255W 


5,600 


2000-11-13 


5,012 


WFPC2 


F336W 


2,000 


2000-11-13 


5,012 


WFPC2 


F439W 


1,200 


2000-11-13 


5,012 


WFPC2 


F555W 


800 


2000-11-13 


5,012 


WFPC2 


F675W 


2,400 


2000-11-13 


5,012 


WFPC2 


F814W 


800 


2003-11-28 


6,122 


ACS HRC 


F625W 


800 


2006-12-06 


7,226 


ACS HRC 


F625W 


1,200 


2009-12-12 


8,328 


WFC3 


F225W 


800 


2009-12-12 


8,328 


WFC3 


F336W 


800 


2009-12-12 


8,328 


WFC3 


F438W 


800 


2009-12-12 


8,328 


WFC3 


F555W 


400 


2009-12-12 


8,328 


WFC3 


F625W 


3,000 


2009-12-12 


8,328 


WFC3 


F814W 


400 



H987 February 23 



Table 2 

HST STIS OBSERVATIONS 



Date 


Days 


sinee explosion^ 


Grating 


Sht width 
(arcscc) 


Slit PA 
(degrees cast of north) 


Exposure time 
(s) 


1999-02-21 




4,381 


G750L 


0.5 




25.6 


10,500 


1999-08-30 




4,571 


G750M^ 


0.1 




27.0 


7,804 


1999-08-30 




4,571 


G750M^ 


0.1 




27.0 


7,804 


1999-08-31 




4,572 


G750M^ 


0.1 




27.0 


7,804 


2004-07-18 




6,355 


G750L 


0.2 




179.7 


5,468 


2004-07-18 




6,355 


G750L 


0.2 




179.7 


5,468 


2004-07-23 




6,360 


G750L 


0.2 




179.7 


5,468 


2010-01-31 




8,378 


G750L 


0.2 




179.7 


14,200 


Note. — The 
and 8,378 are 


slit positions for day 4,571/4,572 are 
shown in Fig. 1111 


shown in 


Fig. [8] and the slit positions 


for days 6,355/6,360 



^987 February 23 

^Central wavelength setting 6581 A 
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